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We have investigated the interplay between the metal-insulator transition and ferromagnetism in 
IIIi-^Mn^V semiconductors. Our study is based on a model in which S = 5/2 Mn local moments are 
exchange-coupled to valence band holes that interact via Coulomb interactions with each other, with 
ionized Mn acceptors, and with the antisite defects present in these materials. We find quasiparticle 
participation ratios consistent with a metal-insulator transition that occurs in the ferromagnetic 
state near x ~ 0.01. By evaluating the distribution of mean exchange fields at Mn moment sites, 
we provide evidence in favor of the applicability of impurity-band magnetic-polaron and hole-fluid 
models on insulating and metallic sides of the phase transition respectively. 

PACS numbers: 75.50.Pp, 75.10.-b, 75.30.Hx 



Introduction- Recent progress [El in the growth of di- 
luted magnetic semiconductors semiconductors that ex- 
hibit ferromagnetism Q at relatively high temperatures 
has suggested exciting new possibilities for devices that 
combine information processing and storage functionali- 
ties in a single material. It is generally accepted that 
Mn acts as an acceptor in these semiconductors, that its 
half- filled d-shell contributes a S = 5/2 local moment 
H-|| to the system's low-energy degrees of freedom, and 
that ferromagnetism is due to interactions between lo- 
cal moments that are mediated by valence-band holes. 
Recent experiments J7|— [Toj] demonstrate that ferromag- 
netism occurs in both metallic and insulating states, and 
that both magnetic and transport properties are sensitive 
to the Mn fraction x and to the density of compensating 
antisite and other defects in the material. The highest 
ferromagnetic critical temperatures (T c ) appear flo|| to 
occur in the most metallic samples. The role of Coulomb 
interactions in the ferromagnetism in these materials is 
subtle. At high carrier densities, well on the metallic 
side of the metal-insulator transition, exchange and cor- 
relation in the hole system is expected to enhance 
ferromagnetism. Well on the insulator side of the tran- 
sition, however, Coulomb interactions create a Mott gap 
that increases the importance of randomness in Mn ion 
positions, suppresses carrier-hopping between Mn sides, 
and eventually turns off the coupling between local mo- 
ments that can lead to ferromagnetism. 

In this Letter we report on a numerical Hartree-Fock 
study of model (III,Mn)V ferromagnets in both metallic 
and insulating regimes. By evaluating the distribution 
of T = Mn exchange mean-fields, defined precisely be- 
low, we find evidence that is generally supportive of the 
impurity-band magnetic-polaron |l2]-|l4| picture that has 
been used to describe ferromagnetism in the insulating 
regime, and of the hole-fluid picture that has been used 
P, pl]p^JT^ | in the metallic regime. Our principle results 
are summarized in Figs. ^ and ^ where we plot vs.T the 



fraction of Mn sites that have exchange mean-fields larger 
than ksT. 
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FIG. 1. f(T), the fraction of Mn that experience a 
mean-field stronger than fcgT, vs. T in the metallic case. 
(Nmti = l.Onm -3 , n h = Nmu/3-) Pnin(n), the probabil- 
ity distribution function for partial hole densities at Mn 
sites in units of the average hole density, is plotted for ma- 
jority (solid line) and minority (dashed line) spins in one 
inset and the quasiparticle density-of-states in the other 
inset. For this case the hole density at an isolated Mn is 
1.63 times larger than the average hole density. The den- 
sity of states D(E) is per occupied hole with energies in 
eV. Note the anomaly at the Fermi energy Ef, indicated 
by a vertical line. The number of disorder realization is 
fifteen for D(E) and ten for f(T) and PMn(n). A a sim- 
ulation cell of side L — 8nm is used. 



In the metallic case 17 (x = 0.05 and p — 3.3 x 



10 cm -3 ) the distribution of coupling strengths is 
peaked at a finite value close to its uniform hole-fluid 
value and all Mn ions are strongly coupled to holes. 
In the insulating case (x = 0.01 and hole density p = 
2.8 x 10 19 cm -3 ), on the other hand, we find that that 
vast majority of Mn spins arc nearly free the few strongly 



coupled Mn ions will form magnetic polarons that grow 
slowly in size and interact more strongly as the tempera- 
ture is lowered. In the following sections we explain how 
these results were obtained and elaborate on their sig- 
nificance for models of diluted magnetic semiconductor 
fcrromagnctism. 
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FIG. 2. As in Figure 1 but for the insulating case. 
(N Mn = 0.25nm~ 3 and n h = N Mn /10). The hole density 
at an isolated Mn site is 19.9 times larger than the average 
hole density in this case. 

Model Hamiltonian.- Most theoretical work on (III,Mn)V 
ferromagnets has started from one of two idealized lim- 
its. Impurity-band models |l2]-|l4|] achieve simplification 
by assuming that the holes that mediate interactions be- 
tween Mn ion moments are strongly localized, whereas 
hole-fluid models achieve simplification by neglecting dis- 
order due to Mn acceptors and other defects and treating 
it perturbatively in estimating transport coefficients. An 
insulator to metal phase transition occurs in these ferro- 
magnets as Mn and hole densities increase; experimen- 
tally |l|] it appears that the metal-insulator transition 
usually occurs near x ~ 0.01, likely depending on the 
density of compensating antisite defects, a quantity that 
is expected to be sensitive the details of sample-growth 
and annealing protocols. Although it is generally agreed 
that an impurity band picture should apply far on the 
insulator side of the transition and a hole-fluid model far 
on the metallic side of the transition, it has not been 
clear which approach is a better starting point in the 
experimentally relevant parameter ranges. We address 
this issue by examining the ferromagnetic ground state 
using a model, and an approximation scheme, that cap- 
tures the physics of both limits correctly. Our study is 
built on k ■ p envelope-function approach description of 
the valence band jL8| and a Hartree-Fock description of 
interactions. 

The four terms in the single-particle part of our band 
electron model Hamiltonian, H° = H K + H Kex - + 
jjMn-h _j_ jjas-h ^ re q mre some discussion. H Kex - repre- 



sents the kinetic exchange interaction between Mn local 
moments, assumed to be aligned in the ferromagnetic sys- 
tem ground state |^9|, H Mn ~ h represents the attractive 
Coulomb interaction between an ionized Mn 2+ acceptor 
and a valence band hole. In an envelope function formal- 
ism, central cell corrections to this interaction are neces- 
sary pel to achieve an accurate description of the isolated 
bound-acceptor limit. H as ~ h describes the repulsive in- 
teraction between holes and antisite defects (represented 
as sites with charge +2). These defects act as deep donors 
partially compensating the Mn acceptors and reducing 
the overall hole density, and provide an important addi- 
tional scattering center. H K is the usual kinetic energy 
term. In this study we ignore |^, [l6| , pl|] mixing between 
heavy and light hole bands by using a simple parabolic 
band approximation. 

+ V( ^U- - V Q e-\ ? -K'?^) + V 2e \ . 

V ef-ifrl ^e\r-R K \ 

(1) 

Here J(r) = ( J pd /(2na 2 ) 3 / 2 ) exp(-r 2 /2a 2 ), r - 
(j x , Ty, t z ) is the Pauli spin matrix vector, I labels Mn 
sites, K labels antisite defect sites, and f2j is the orien- 
tation of the I th Mn spin. The term in the potential 
proportional to Vo is the central cell correction. Both 
Mn ions and antisites were distributed randomly |2^] 
in a cube of side L. The long range of the Coulomb 
interaction requires overall charge neutrality so that 
nh — Nmu + 2N as — 0, where nh is the density of holes, 
Nmu the density of Mn ions and N as the density of an- 
tisites. In the ground state, all Mn spins are oriented 
along the z direction and H° is block diagonal in its spin 
indices. In the spin wave configurations discussed below, 
x and y components of the spins are present, doubling 
the dimension of the matrix that must be diagonalized 
numerically. The wavevector cutoff in the quasiparticle 
wavefunction expansion [ p3| was tested by computing the 
binding energy of a hole to an isolated Mn, comparing 
with the results of Bhattacharjee and Benoit a la Guil- 
laume J2(J who find that a binding energy of 112meV, 
86 meV when the kinetic exchange term is neglected, 
and 68meV when the central cell correction is also ne- 
glected. Our results are 124meV, 88 meV, and 48meV 
respectively, demonstrating an adequate description of 
the completely isolated Mn limit. 

In order to capture the correct physics of both metal- 
lic and impurity band limits, hole-hole interactions must 
be described using an approximation that accounts for 
screening effects in the metallic regime and avoids artifi- 
cial self-interaction effects in the localized regime, moti- 
vating our use of Hartree-Fock (HF) theory |24| . The HF 
quasiparticle Hamiltonian is H HF = H° + (V H + V x ) 
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where the Hartree and Fock matrix elements, V- - and 

k(7,k f cr' 

g , can be expressed in terms of the density ma- 
trix. We evaluate the density matrix in the fc-space rep- 
resentation used to expand our envelope function quasi- 
particle wavefunctions: pr n , = Y"! n«c- c- *, where 

| a) = J^kcr C ^\k a )' an d n a is a quasiparticle occupation 
number. Our HF scheme becomes exact in the strongly 
localized limit, since a localized quasiparticle does not 
interact with itself, but is expected to overestimate inho- 
mogeneity in the metallic limit because it neglects quan- 
tum fluctuations in the many-hole state. 
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FIG. 3. Participation ratios vs. energy in a simulation 
cell of side L = 8nm for majority (solid line) and minor- 
ity (dashed line) quasiparticles. The Fermi energies are 
indicated by solid vertical lines. The upper panel is for 
the metallic case of Fig. 2 and the middle panel is for the 
insulating case of Fig. 3 and these results were obtained 
by averaging over 15 disorder realizations. The bottom 
panel shows the system size dependence of majority-spin 
participation ratios at the Fermi energy for the two cases. 

Participation ratios.- In order to verify that our model 
correctly describes the metal-insulator transition in these 
materials, we measured the localization properties of our 
Hartree-Fock quasiparticles by evaluating their partici- 
pation ratios, P a = 1/(L 3 J d 3 r|* Q (r)| 4 ), where * a (r) = 
(r\a) is a normalized quasiparticle wavefunction. (For a 
localized state P a ~ (£ a /L) 3 , where £ Q is the localization 
length.) Fig. 3 compares participation ratios evaluated 
at x = 0.05 and 0.0125, assuming compensated carrier 
densities = Nmu/3 and nh — Nmu/IO, respectively. 
The participation fractions extrapolated to infinite vol- 
ume are clearly finite and therefore consistent with metal- 
lic transport at x — 0.05 whereas they are consistent with 
zero and localized quasiparticles at x = 0.0125, in agree- 
ment with experiment. 

Impurity- Band Picture vs. Hole-Fluid Picture.- The ex- 
change mean-fields used to construct Fig. [l] and Fig. || 
are given by H effJ = J drJ(r- #j)(nj(r) - n T (r))/2 w 

Jpd(n>i(Ri) — n ~\(Ri))/2 where the partial densities are 



determined by solving the Hartree-Fock equations self- 
consistently. In the extreme impurity band limit, H e ff j 
would have the value J p d n max/2 ~ 25meV for a frac- 
tion rih/NjMn of the local moments and much smaller 
values for all others. Our findings are in qualitative 
agreement, although the peak H e ff value is somewhat 
smaller and there is no sharp peak in the distribution 
function near the peak value, presumably because of 
the variable Mn and antisite defect environment ex- 
perienced by localized holes. (The shallow impurity 
Bohr radius a* B calculated using the heavy-hole mass 
is ~ l.Onm; the peak hole density, 4.9 x 10 20 cm~ 3 is 
slightly larger than the shallow impurity value because 
of central cell corrections.) Kaminski and Das Sarma 
|fl3f have recently estimated the impurity-band mag- 
netic polaron model T c , relating it to Nmu, nh and 
the maximum mean-field exchange coupling by T c ~ 

0.5(WWn fc ) 1 /*(n hO l» 3 ) 1 /e ff(na!e exp(-0-W/nh 1/3 a%)/k B 
Naively applying this formula to our x = 0.01 case with 

the shallow impurity Bohr radius yields T c ~ 15A'. Note 
that T c is expected to decrease rapidly at lower hole 
densities because of the weakening magnetic polaron 
coupling. These estimated T c 's are somewhat larger 
than what has typically been observed at x ~ 0.01 in 
(Ga,Mn)As, possibly because the theory neglects Pauli 
exclusion principle effects, superexchange interactions, 
that favor opposite orientations of nearby impurity-band 
holes and oppose ferromagnetism. Nevertheless, it ap- 
pears that this picture provides an excellent qualitative 
description of ferromagnetism at x ~ 0.01 and smaller. 

In the metallic case, the exchange mean-field distri- 
bution differs qualitatively, with strong couplings com- 
mon and few weakly coupled moments. Our results 
likely overestimate the degree of density variation in sam- 
ples with x ~ 0.05, both because we have neglected 
Mn acceptor-antisite correlations (22) and because of our 
use of the Hartree- Fock approximation. (The corre- 
sponding self-consistent Hartree calculations lead to more 
sharply peaked Mn site majority-spin density distribu- 
tions.) We note that the most likely hole density at a 
Mn site is approximately 1.5 times larger than the av- 
erage hole density and that the most likely mean-field 
coupling strengths are again somewhat larger than the 
uniform hole-fluid value Hfi = nJ p d/2 w 190K. 

There are important differences between metallic and 
insulating cases in the physics that controls T c . The 
metallic T c can be limited either |2^,|2^] by the system's 
stiffness against collective magnetization orientation vari- 
ation, or if the stiffness is large enough by the compe- 
tition between local moment entropy, exchange interac- 
tions, and the band energy cost of hole-spin polarization 
that is captured by the hole-fluid mean-field-theory T c 
expression p6fl . We have estimated the magnetic stiff- 
ness of both insulating and metallic states by evaluating 
the HF electronic energy cost of an imposed spin-wave 
with wavevector Q. If disorder were neglected, the pro- 
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cedure we follow in evaluating spin spin-wave energies 
differs from the theory |^7j of Konig et al. only through 
retardation effects. 
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FIG. 4. Average spin-flip energy vs. wavevector. The 
upper curve is for the metallic case of Fig.l while the lower 
curve is for the insulating case of Fig. 2. The wavevector 
is normalized to a Debye wavevector defined by the Mn 
density in each case, k c = (6'7r 2 iVjvi n ) 1/ ' 3 . The inset shows 
a magnified view of the small wavevector regime. 



Figure 4 displays average spin reversal energies for metal- 
lic and insulating states. Note that the energy cost of 
slow (long wavelength) spin-direction variations is ex- 
tremely small in the insulating case because these are 
sensitive mainly to the weak interactions between mag- 
netic polarons that are necessary for long-range magnetic 
order. It follows that mean-field will completely fail in 
estimating the critical temperature of insulating DMS 
ferromagnets. The energy cost for spin-direction varia- 
tions that we find in the metallic case is in qualitative 
agreement with hole- fluid model results, although it is 
larger at short wavelengths. In the hole-fluid mean- field 
theory the energy cost of spin-direction variations would 
equal Hfi at all wavevectors. It follows that for the 
parabolic band model we have studied substantial cor- 
rections p8| , p5[ would apply to its mean-field T c . It is 
important to note, however, that valence band spin-orbit 
interactions stiffen collective spin-variations in these fer- 
romagnets p£| , reducing corrections to mean-field-theory 
T c estimates. 
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